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Abstract. We investigate the derealization transition appearing in an exclusion process with two internal 
states resp. on two parallel lanes. At the transition, delocalized domain walls form in the density profiles of 
both internal states, in agreement with a mean-field approach. Remarkably, the topology of the system's 
phase diagram allows for the derealization of a (localized) domain wall when approaching the transition. 
We quantify the domain wall's derealization close to the transition by analytic results obtained within the 
framework of the domain wall picture. Power-law dependences of the domain wall width on the distance 
to the derealization transition as well as on the system size are uncovered, they agree with numerical 
results. 

PACS. 5.40.-a, 05.60.-k, 64.60.-i, 72.25.-b 



1 Introduction 



Driven one-dimensional transport phenomena pQ currently 
receive much attention as they constitute a challenging 
class of non-equilibrium dynamical systems. In such sys- 
tems, collective effects induce unexpected phenomena, in- 
cluding e.g. boundary-induced phase transitions [5] or pat- 
tern formation [3]. Possible applications are found in a 
large variety of contexts, ranging from biology (e.g. the 
motion of ribosomes along mRNA [Jj or molecular mo- 
tors on intracellular filaments 5,6, 7,8, 9, 10 , 11) to elec- 
tron hopping transport with applied voltage |12|J13j and 
■ vehicular highway traffic [14,15]. 

As unifying descriptions of such non-equilibrium sys- 
tems are still lacking, much effort is devoted to the under- 
standing of particular models and the identification of uni- 
versal properties as well as analytic methods. In this con- 
text, the Totally Asymmetric Exclusion Process (TASEP) 
has emerged as a paradigm (see e.g. [T5UT7] for a review). 
There, particles move unidirectionally from left to right 
through a one-dimensional lattice. The injection resp. ex- 
traction rates at the left resp. right boundaries serve as 
control parameters; tuning them, different phases of the 
stationary-state density are observed. Although being ex- 
actly solvable [T8l[T9l[20] . a mean-field (MF) approach al- 
ready yields the exact phase diagram, originating in an 
exact MF current-density relation [TB]. Beyond MF, fluc- 
tuations have been successfully taken into account in the 
framework of the domain wall picture |21,22J. There, one 
considers a domain wall separating a low-density (LD) re- 
gion from a high-density (HD) one. The domain wall's dy- 



namics yields information about the phase behavior and 
boundary effects arising in finite systems. 



Fig. 1. (Color online) Illustration of an exclusion model with 
two internal states, denoted as spin-up (j) and spin-down (J.). 
Particles are injected at the left boundary at rates a\a^ , and 
extracted at the right at rates /3 T , (3^ . Within bulk, particles 
move unidirectionally to the right, or flip their state (at rate 
oj), always obeying Pauli's exclusion principle. 



Recently, we have proposed a generalization of TASEP 
where particles possess two weakly coupled internal states 
[T3"ll23j . As an example, these internal states may corre- 
spond to different parallel lanes in vehicular traffic on 
highways [14] . Concerning intracellular transport, micro- 
tubules consists of typically 12-14 parallel lanes, and molec- 
ular motors progressing on them may (though rarely) switch 
between these lanes. Our work provides a minimal model 
that takes such lane changes into account. Also, the inter- 
nal states may describe spin states of electrons, e.g., when 
performing hopping transport through a chain of quantum 
dots [12] , suggesting possible application of our model to 
spintronics devices. 

In the subsequent paper, we adopt the language of 
spins for the internal states. The system's dynamics, de- 
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scribed in the following, is depicted in Fig. [T] Particles 
with spin-up (down) state enter at the left boundary at 
rate (a^), under the constraint of Pauli's exclusion 
principle. The latter means that each lattice site may at 
most be occupied by one particle of a given state, such that 
particles with spin-up and spin-down may share one site, 
but no two spin-up or spin-down particles are allowed. (In 
the context of multi-lane traffic, this translates into sim- 
ple site exclusion.) Within the lattice, particles hop to the 
neighboring right lattice site at constant rate (which we 
set to unity), again respecting Pauli's exclusion principle. 
Spin flip events may occur, we denote the respective rate 
by uj. Finally, particles are extracted at the right boundary 
at rate ft (P^), depending on their spin state. 

In p~3l[23] we have taken advantage of a MF approach 
in a continuum limit, see e.g. [24] for a review, to de- 
scribe the emerging non-equilibrium stationary density 
profiles. Comparison to data from stochastic simulations 
has, through finite size scaling, uncovered an apparent ex- 
actness of the analytic results in the limit of large sys- 
tem sizes. We have traced back this exactness to the weak 
coupling of the two internal states, see below, as well as 
the exactness of the MF approach for TASEP [IB]. As it 
implies the exactness of the analytically derived phase di- 
agrams (see Fig. [2] for a two-dimensional cut), we have 
full knowledge of the system's phase behavior. The latter 
exhibits a rich variety of phases, in particular, a localized 
domain wall, separating a low-density (LD) from a high- 
density (HD) region, may emerge in the density profile of 
one spin state. We refer to this situation as a coexistence 
of LD and HD, abbreviated as LD-HD. Continuous as well 
as discontinuous transitions between the different phases 
appear and induce multicritical points (e.g. .Ain and .4ex 
in Fig. Hereby, as we already announced in [T3"l[2"3"] . 
domain wall derealization appears at the discontinuous 
transition. Namely, as an example, consider the red (grey) 
line in Fig. [2j It crosses a discontinuous transition, which 
seperates two LD-HD phases. In both, localized domain 
walls emerge, but at different positions in bulk. Upon 
crossing the derealization transition along the red line, 
the domain wall delocalizes, and then localizes again at 
a different position. The aim of this article is the inves- 
tigation of the domain wall's derealization process when 
approaching the discontinuous transition. 

In this paper, we present detailed investigations on 
the derealization transition. Starting from MF consid- 
erations, we show that at the transition line delocalized 
domain walls appear in the density profiles of both spin 
states. More precise, the MF approach predicts a one- 
parameter family of solutions with domain walls for spin- 
up and spin-down, which we conclude to perform coupled 
random walks. Further on, including fluctuations via the 
domain wall picture allows to quantify the derealization 
of a localized domain wall when approaching the dereal- 
ization transition, e.g. along the path depicted in red (or 
grey) in Fig. O 

The outline of this paper is the following. In the next 
section we review the MF approach (see our previous ar- 
ticle [23 for a detailed discussion), and present the solu- 
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Fig. 2. Phase diagrams for the density of spin-up (a) and spin- 
down (b): two-dimensional cuts for fixed values of = 0.25, 
ft — 0.3 and in the asymptotic limit of large Q. Hereby, 
O — uj/L is the gross spin flip rate. Lines of continuous tran- 
sitions (thin) intersect the derealization transition line (bold) 
in multicritical points ,4in and *4ex- The scope of this paper is 
to elucidate the system's behavior at the derealization tran- 
sition as well as when approaching the latter along the path 
shown in red (or grey). 



tion for the density profiles at the derealization transi- 
tion line. The solution turns out to be not unique, but 
is a one-parameter family, indicating the derealization. 
Our stochastic simulations are described in Section [3] Sec- 
tion 3] goes beyond MF by taking fluctuations into ac- 
count. Within the domain wall picture, we investigate the 
derealization of a domain wall when approaching the de- 
localization transition. As a result, we find power law de- 
pendences for the width of the domain wall distribution 
on the system size and the distance to the derealization 
transition line. Our conclusions are presented in Sec. [5] 



2 Predictions from the mean-field approach 

We are interested in the stationary density profiles which 
emerge in the system's non-equilibrium steady state. There- 
fore, denote the occupation number of site i for spin-up 

resp. spin-down state by nj resp. n\, i.e. G {0, 1}, de- 
pending on whether this site is occupied by a particle with 
corresponding spin state or not. Performing sample aver- 
ages, we obtain the average occupation, pj = (nj ). In 
the mean-field (MF) approximation, higher order correla- 
tions between the occupation numbers are neglected, i.e. 
we assume 



(i) 



The dynamical rules of the system lead to equations for 
the densities in the stationary state; see Ref. [23] for a 
detailed discussion. Further on, in a continuum limit, the 
total length of the lattice is set to unity and the limit 
of number of lattice sites L — > oo is considered. Hereby, 

the densities p| approximate smooth functions p^^(x) 
with x G [0,1]; for the latter, we eventually obtain two 
coupled differential equations: 



cV T - n{p l - p ] ] , d x j l = rV - P 



(2) 
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Here, we have defined currents for the individual spin 
states: j^^(x) — p^^(x)[l — p^^(x)]. Also, the gross 
spin flip rate J? = luL was introduced; in a mesoscopic 
scaling, we keep this rate fixed when performing the limit 
L — > oo, see Refs. 25,26 . In this way, competition be- 
tween the bulk processes (spin flips) and the boundary 
processes (particle injection and extraction) emerges. This 
weak coupling is opposed to strong coupling, where u> is 
kept constant when considering large systems, which leads 
to different effects [2"71l2"51l2U] . A closed analytic form for 
the solution to Eqs. © is feasible and has been presented 
in [53]. 

The injection and extraction processes lead to bound- 
ary conditions for Eqs. @: 



P l (0) = 



T 

err ' 



*eff> 



off ' 



(3) 
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Fig. 3. Domain walls in the densities of both spin states at the 
derealization transition: results from the MF approach. The 
positions of the domain walls are only determined up to one 
degree of freedom: choosing the position zj, for the domain wall 
of spin-up yields the position xj, for spin-down and vice versa. 
The densities found in stochastic simulations are the average 
over this one-parameter family of solutions. Parameters are 
a T = 0.45, a 1 = 0.2, /3 T = 0.36, /3 l = 0.23 and Q = 0.3. 



where we have introduced effective rates a 
cording to 



ra- 
ce . 



old) 

' PcS 



(4) 



They reflect the fact that bulk processes limit the individ- 
ual spin currents to maximal values of 1/4, corresponding 
to densities of 1/2. Injection or extraction rates exceeding 
this value cannot lead to larger currents, but effectively 
act as 1/2 (see [23]). 

The boundary conditions ([3J apparently overdetermine 
the system of two differential equations ([2]). Indeed, two 
types of singularities may occur where the densities ex- 
hibit discontinuities: boundary layers (discontinuity at the 
boundary) and domain walls (discontinuity in bulk). An 
analytic description of the domain wall positions is feasi- 
ble from the observation that the spin currents j^^(x) 



situation, delocalized domain walls appear in the density 
profiles of both spin states. Indeed, only when Jin = Jex, 
density profiles satisfying the boundary conditions at the 
left boundary as well as at the right one are feasible. In- 
stead of discontinuities at the boundaries, they exhibit 
discontinuities within bulk: domain walls. In the follow- 
ing, we want to present the MF picture for the density 
behavior at the dclocalization transition and show that a 
one-parameter family of solutions emerges within the MF 
approach. 



2.1 The one-parameter family of analytic solutions 

The generic picture of the MF analytic solution at the de- 
localization transition is presented in Fig. [3J As the equa- 
tion Jin = Jex is fulfilled, the density profiles can match 
the left as well as the right boundary conditions. This 
implies that, in bulk, domain walls arise at positions x\ 



(x)[l 



P 



ru) 



(x)] at these positions must be continu- and x l in spin-up resp. spin-down state; we have chosen 



ous. The only solution to this condition is that a density 
pTU) changes to the new value 1 — pTU) a t the position of 
the discontinuity. 

As particles cannot leave the system in bulk, the parti- 
cle current J = p +j± is spatially conserved. In case that 
no boundary layer occurs at the left, i.e. the densities do 
not have any discontinuities there, it is given from (J3j by 

the value Jin = a lff{^ — a lff) + a cff(l — a c«)- ^ n this case, 
the system is dominated by injection processes; we refer 
to the region in parameter space, where this behavior oc- 
curs, as the IN-region. On the other hand, also extraction 
can determine the system's behavior, the particle current 
is then given by J EX = /3 e T ff (l - /3 T ff ) + #^(1 - /3 c X ff ), occur- 
ring in the EX-region. 

The delocalization transition occurs at Jin = Jex- 
There, the system is in a superposition of the injection 
dominated and the extraction dominated behavior. In this 



< x+j, the other case is obtained by the symmetry 
of the two spin states. In the vicinity of the left bound- 
ary, the densities are given by the analytic solution p^ 
obeying the left boundary conditions pj(x = 0) = al s , 
pj(x — 0) = oigjj. At the point x^, the density of spin-up 
jumps from the value p[{x\ u ) to the value 1 — /^(x^), while 
the density of spin-down is continuous. To the right of x^, 
the densities follow a different branch of analytic solution, 
we name it p]^ as the intermediate branch. It is deter- 
mined by the boundary conditions p^x^) = 1 — pl(x1,) 
and Pmix)^) = p\{xl w ). At the point x^, the density of 
spin-down changes discontinuously from a value p\ n {x\ J) ) 
to a value 1—p^x*,), while the density of spin- up remains 
continuous. To the right of x^, the densities obey the an- 
alytic solution pl^ satisfying the boundary conditions at 
the right: pi {x = 1) = 1 - /3 c T ff , p\.{x = 1) = 1 - ^s- 
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Fig. 4. (Color online) Analytic versus numerical results for the 
density profiles at the delocalization transition. The analytic 
approach predicts a one-parameter family of solutions, we show 
the one with the minimal domain wall positions (red or grey) 
and with the maximal one (blue or dark grey) as well as an 
intermediate (light blue or light grey). Stochastic simulations 
yield an average over these solutions, resulting in densities in- 
terpolating between the extrema (black). Parameters are the 
same as in Fig. [3] 



Now, for the situation to be feasible, we have additional 
conditions coming from considering the continuity of the 
spin currents at the point x^,, namely pl(x*,) — /0„(x^) 
and p*(x^) = 1 — Pmi x w)- As the particle current is con- 
served, these two conditions are not independent: the va- 
lidity of one implies the validity of the other. Only one 
condition thus exists for two variables, the latter being 
the positions of the domain walls xj, and x^ . Our consid- 
erations therefore leave us with a one-parameter family of 
solutions; as parameter, we can e.g. take one of the domain 
wall positions, x]^ or x*,. 



The boundaries for the residence regions of the domain 
walls originate in the constraint that the domain wall po- 
sitions must lie between and 1. Indeed, let us consider 

the minimal positions xl} mm ^ , xi which are feasible 
for the domain walls in the spin-up and spin-down state. 
We had found in Subsec. [2] that fixing the domain wall 
position in the density profile of one of the states deter- 
mines the position of the other domain wall. The situation 



i(min) 



™T(min) 
w 

„T(min) 



= can thus only emerge if x 



i(min) 



= 



X LI 

induces Xw' ULU> = (corresponding to the multicritical 
point Am, see the detailed discussion in Ref. [23]). In gen- 
eral, only one of these minimal domain wall positions is 
at 0, and induces a minimal position for the other one 
which is larger than 0. Analogously, in general, a maximal 
domain wall position of 1 only occurs for one of the spin 
states, and induces a domain wall position smaller than 1 

for the other. The exceptional case xi| max ' ) = xlj m = 1 
occurs at the multicritical point -4ex, as the latter lies at 
the transition from a coexistence phase to a HD phase, for 
both spin states, see Ref. [23] 

Which situation occurs, may be read off from the phase 
diagrams by considering the phases in the vicinity of the 
delocalization transition. As an example, we consider a 
point on the delocalization transition line in Fig. [2] (a), 
(b), at the intersection with the path shown in red. In 
its vicinity, in the EX-region, the density of spin-down 
is in a pure HD-phase, such that a domain wall at posi- 
tion x!^ = forms there when approaching the delocal- 
ization transition. The density of spin-up exhibits a local- 



ized shock at a position < x 
xi, (min) = and < xl {mia) < 
the IN-region in the vicinity of this point yields the maxi 



, < 1. We thus observe 
1. Similarly, considering 



J. (max) 



= 1 and 



mal feasible domain wall positions there: 

o < xi (max) < i. 

The analytic solutions to the density profiles in the sit- 
uation of xl} mm ^ , Xw™'^ (red) as well as for xjj m , Xw™^ 
(blue) are shown in Fig.[U The intermediate solution (light 
blue) varies between these two extrema. In the system, av- 
eraged over sufficiently long time, this results in smoothened 
density profiles. In Fig.[¥]we show the results from stochas- 
tic simulations as black lines, agreeing with our analytic 
considerations. 



3 Stochastic simulation methods 



2.2 Spatial boundaries for the domain walls 

In the course of time, the system takes all states specified 
by the one-parameter family of solutions. Averaged over 
sufficiently long times, the densities represent the average 
over the latter, resulting in a smoothening of the density 
profiles. In Fig. [4] we show how the results may look like. 

We observe that certain spatial boundaries (x^""' and 

it"™' in the above picture) exist. They define the bulk 
region where the domain walls may appear. The descrip- 
tion of these boundaries is the scope of this subsection. 



To validate our analytic calculations, we have carried out 
extensive stochastic simulations. An efficient simulation 
method originally due to Gillespie [30][3"T] was implemented. 
There, in each time step, a random number determines 
whether particle injection, particle extraction, a spin flip 
event or a particle hopping forward may occur in the next 
time step. The time interval to the next process is chosen 
from an exponential waiting time distribution. 

Simulations in the neighborhood of the delocalization 
transition need long waiting times until the densities reach 
stationary profiles, due to the random walks performed by 
the domain walls. We carried out simulations with up to 
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Fig. 5. The width a of the domain wall distribution in the 
spin-up density profile for the case of a localized domain wall 
there, at a small distance S to the derealization transition. 
Numerical results are shown in a double logarithmic scale, de- 
pending on |Ld^|. The data are obtained from system sizes of 
L = 5000 (open circles) and L = 10000 (open boxes), with d)~ 
varied between 1 x 10 -4 and 1 x 10 -3 . The solid line indicates 
the slope — |. Approximately, we recover a ~ \LS\^ 1 ^ 2 . 



10 s time windows, each consisting of 10 x L steps of up- 
dating. An example of resulting density profiles is given in 
Fig-Sfor a point on the derealization transition line, the 
densities are observed to interpolate between two extremal 
analytic solutions. 

In the subsequent section, we investigate the scaling 
of the width a of the domain wall when approaching the 
derealization transition. For the numerical data, we have 
computed the stationary density profiles for system sizes 
of L = 5000 and L = 1000 at different distances to the 
derealization transition. The width a was obtained by 
fitting the densities profiles in the vicinity of the domain 
wall with the function A ■ erf[(x — b)/a] + C + Dx, with 
parameters A, b, C, D, a. The first term describes a Gaus- 
sian distribution of the domain wall, while the latter two 
terms account for the background density profile (to lin- 
ear order in x) . The results are shown in Fig. O and agree 
well with the analytic predictions. We attribute deviations 
to the fact that for large a the above fitting becomes less 
accurate as the domain wall distribution deviates from a 
gaussian. 



4 Beyond mean-field: Approaching the 
derealization transition 



In the previous Sections, we have described how delocal- 
ized domain walls in the density profiles of both states ap- 
pear at the derealization transition. Here, we investigate 
the emergence of derealization when approaching the de- 
localization line by applying the domain wall picture. In- 
cluding fluctuations, the latter allows us to calculate the 
width of the domain wall distribution for finite system 
sizes L. It diverges as one approaches the derealization 
line, yielding delocalized domain walls. 



4.1 The domain wall picture 

The domain wall picture [21] allows to include fluctu- 
ations, such as the particle- number fluctuations in the 
TASEP [22], helping in understanding its dynamics [32], 
and thus go beyond the MF approximation. In this ap- 
proach, we start with the assumption that sharp and lo- 
calized domain walls appear in the density profiles of both 
states. Injection and extraction of particles as well as spin- 
flips induce dynamics for them, they perform random walks. 
Our focus is on the emerging stationary state. As the lat- 
ter arises in a non-equilibrium system, it need not obey 
detailed balance. Writing down the master equation for 
the random walk, we are able to calculate the fluctuations 
of the domain wall positions around the stationary values, 
yielding the width of the domain walls or boundary layers 
which form in finite systems. 

The dynamics induced by entering and exiting parti- 
cles or spin-flips is the following. If the total number of 
particles of given internal state in the system is increased 
by one, the domain wall in the density profile of this state 
moves a certain value to the left. In the other case, if the 
number of particles of given state is decreased by one, the 
domain wall of this state moves to the right. 

Denote A^ resp. A± the height of the domain wall in 
the density profile of the spin-up resp. spin-down states. 
They are functions of the domain wall positions: A^ = 
A^ (xl„, x*,), Ai = Ai(xlj,x*,). Increasing the number of 
spin-up particles by one shifts the position x]^ by the value 



5x1 



-L ZL , and decreasing its number results in a 



shift Sx^ = L 1 A^ 1 . The changes Sx^ result analogously 
from particles with spin-down entering and exiting. 



4.2 Fokker-Planck equation and the generic form of 
the domain wall distribution 

In our model, six processes alter the particle number of 
one or both internal states: 

(i) A particle with spin-up entering at the left bound- 
ary 

(ii) A particle with spin-up leaving at the right 
boundary 

(iii) A particle with spin-down entering at the left 
boundary 

(iv) A particle with spin-down leaving at the right 
boundary 

(v) A particle with spin-up flipping to spin-down 
in bulk 

(vi) A particle with spin-down flipping to spin-up 
in bulk 

Their probabilities and resulting changes Sxj^, Sx^ are 
listed in Tab. Q] 

To shorten our expressions, we introduce the vectors 
x w = 04,,a4) and Sx w = (8x1, Sx^). Defining P{x w ,t) 
as the probability density of finding the domain wall po- 
sitions x]^, X+, at time t, the processes (i)-(vi) allow us to 
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Process 


Probability W 


(5r T 


(St 1 


(0 


°ia(l ~ Q cff) 







(ii) 









(iii) 









(iv) 


/3i ff (i - /&) 







(v) 


^/ i P T (« ; )[l-P i (a ; )]da: 






(vi) 


rtJ- n V(*)[l-P T (*)]<*= 







Table 1. Six processes lead to changes in the domain wall posi- 
tions. This table lists their probabilities and resulting changes 
6x w , Sx w . 



write down the master equation 
d t P{x w , t)=Y^ { p (xw + s %w, t)W(x w + Sx w -> x w ) 

(5) 

x x + Sx w ) are 



5x„, 



The transition probabilities W(x u 
the ones listed in Tab. [1] 

For the expectation values of the changes Sx^ , 5x*, we 
obtain 

(5x1) =i-M- 1 {/3 o T ff (i _ fa) _ a T ff(1 _ a T ff) 

+ n f [p\x) - P l (x)]dx} , 

Jo 

(Sxl) =i-M- 1 {^ ff (l - /4) - ai s (l - ai s ) 

+ Q I [p l (x) -p\x)]dx} , (6) 
Jo 

which, of course, depend on the densities p^ and p- 1 and 
thereby on the domain wall positions x^, x*,. Investigat- 
ing the fixed point of the random walk, i.e. the values 
xJu, x^ where (Sxj^) = (5x^) — 0, we find as a necessary 
condition al s (l - afc) + 0^(1 - a l eS ) = /3 c T ff (l - /3 e T ff ) + 

/3g ff (l — /?off)i or ^in = Jex, describing the delocalization 
transition. Thus, domain walls within the density profiles 
of both states are only feasible there; otherwise, at most 
one of the domain wall positions, say x]^, can lie inside the 
bulk. The other one, say x^, is driven outside of the sys- 
tem and turns into a boundary layer. We then set x^ = 
or = 1, depending on whether a HD or a LD phase 
occurs. With this convention, the values x]^, x^ are given 
by the MF analytic solution. 

We are interested in the fluctuations of the domain wall 
positions around the mean values x^, x^, and therefore 
define the quantities 



T T 

y = x w 



w ' 

1 



(7) 



as the deviations. Again, they are arranged in a vector 
y = (y,y^), and we have Sy = Sx w . 



Applying the Kramers-Moyal expansion [33j of the mas- 
ter equation (O around the mean-field values to second 
order in the quantities y^ , results in the Fokker-Planck 
equation 

d t P(y, t) = -d^(y)P(y, *)] + -didj [Bij (y)P% t)] . (8) 

Here, the indices i,j stand for spin-up and spin-down. 
In the above equation, the summation convention implies 
summation over them. The partial derivative di is the 
short-hand notation of d/dy 1 . 

The quantities a.i and By are, according to the Kramers- 
Moyal expansion: 



<k{y) 



5y l W(y 



B^(y)=J2 6 y ls y 3w (y 

Sy 



y + Sy), 
->y + 5y). 



(9) 



For our case, they are given in App. [XJ 

Now, we expand the Fokker-Planck equation around 
y = 0, using the linear noise approximation of van Kam- 
pens ^-expansion [33]. As B^(y = 0) ^ 0, we consider 
approximately Bij(y) Bij(y = 0) = By. On the other 
hand, aj(y = 0) = is possible, and we include the first 
order: 

ai(y) « a,(y = 0) + y 3 djai(y)\ . (10) 



Defining A io = d 3 a t (y)\ y=Q and a, ; 
Fokker-Planck equation |8]) turns into 



'Mr 
a>i(y 



d t P(y, t)=- a%P(y, t) - [A ij y :j P(y, t)] 



-B^P^t) 



0), the 



(11) 



The above expanded Fokker-Planck equation may be 
solved by the ansatz 



p(y, t) ~ exp [ - ±(£-%-2/y + erV] 



(12) 



which is the generic form of the domain wall resp. bound- 
ary layer distribution. When put into Eq. fill) , it leads to 
a matrix equation for S: 

SA T + AS + B = 0, (13) 

the solution is given in App. [X] £ is then determined by 

tr^-iS-^Ajfa. (14) 

Here, we continue by discussing the physical meaning 
of S and £. First, assume that a domain wall emerges 
in the density profile of the spin-up states. Then, x^ lies 
in bulk, and we have a-f = 0. According to Eq. (fl"4"]) . this 
implies f;^ 1 = 0. The domain wall distribution in the spin- 
up density is therefore Gaussian and its width given by 

(S TT )V*. 

Second, for a boundary layer forming in the density of 
spin-up states, it follows that a-f 7^ and therefore also 
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£| 1 ^ 0. As for large enough systems fluctuations are 
small (y , Jr <C 1), the contribution to the domain wall 
distribution coming from ^^y^ is the dominating one. 
The value £j thus describes how far the boundary layer 
extends into bulk, we refer to it as the localization length. 



4.3 Approaching the coexistence line 

When the derealization line is approached, we know from 
the MF analysis that domain walls delocalize. In the do- 
main wall picture, this must result in a diverging width 
of the domain wall resp. boundary layer distributions. In 
this section, we want to show that this divergence indeed 
emerges within the above description. Also, we calculate 
the corresponding exponents and compare our findings to 
stochastic simulations. 

We define the quantities 

=/3 c T ff (l - /3 c T ff ) - aUl - al s ) + Q / V - P l ]dx , 

Jo 

d L =/4(l - /4) - a c V 1 - a l cS ) + n f \pl- p^]dx , 

Jo 

(15) 

as a measure of distance to the derealization transition 
line, where S = — is encountered. Our aim is to 
calculate S and £ for d',d^ 0, i.e. along a path that 
ends at the coexistence line. 

Different paths are possible. Applying the two symme- 
tries of the model, they fall into two classes. First, we may 
have a localized domain wall in the density of one of the 
internal states, and a pure LD or HD phase for the other 
one. In the second class, the densities of both states are 
in pure LD or HD phases. 

Note that the second case is similar to the situation in 
TASEP. There as well we may approach the coexistence 
line on a path along which LD or HD is encountered; in 
this situation, the localization length £ diverges. In our 
model, exactly the same behavior emerges: the localiza- 
tion lengths describing the boundary layers in the 
densities of spin-up resp. spin-down state, both diverge. 

The first case, corresponding to the path shown in red 
in Fig. [21 does not possess an analogy to earlier studied 
ASEP models. Also off the derealization transition line, 
we have a domain wall, which is localized. It allows us to 
study how its width diverges when approaching the IN- 
EX-boundary, i.e. how its derealization arises. 

We start our considerations with the first case. As a 
representative example, a situation with a localized do- 
main wall in the spin-up density and a pure LD or HD 
phase for spin-down is studied. This case implies d' = 
and di = /3 c T ff (l - Pl s ) + (3 l cS (l - (3 l cS ) - - c^) - 

cVg ff (l — Qfg ff ) ^ 0. Both rfJ- < and d 1 - > are possible. 
Our focus is on the width a — (Zfj) 1 / 2 of the domain 
wall distribution in the spin-up density. More specific, we 
aim at finding the scaling behavior of a depending on the 



system size L and the distance d^ . We expect that a — > 
for increasing system size L — ► oo on the one hand, and 
that a diverges for d+ — > on the other hand. 

In App.[Aj we solve the expanded Fokker-Planck equa- 
tion (jlip . Considering the solution (|2l))) . we recognize that 
£ - (det^)- 1 , and from ([22]) we infer det A - d l . Com- 
bining these two results, we arrive at a = (Zjj) 1 / 2 ~ 
jrfJ- 1 -1 / 2 , such that a diverges when the derealization tran- 
sition line is approached. The mathematical reason is that 
one eigenvalue of A goes to zero (such that det A — > 0). Of 
course, the eigenvector to the vanishing eigenvalue corre- 
sponds to the direction along the one-parameter curve on 
which the mean-field values xl^ , can reside when the 
derealization transition line is reached. Thus, the diver- 
gence of the width of the domain wall distribution orig- 
inates in the one degree of freedom which exists at the 
derealization transition. 

The dependence of a on the system size L is also of in- 
terest. Remembering A ~ L^ 1 and B ~ L~ 2 , and with 
help of Eq. ([26| . we arrive at a ~ L^ 1 / 2 . Thus, with 
increasing L, the domain wall distribution gets sharper 
aligned to its mean value. Together with the previous re- 
sult, we obtain the following scaling behavior for the width 
a = (Zji) 1 / 2 of the domain wall distribution in the spin- 
up density: 

a ~ \Ld l \~ 1/2 . (16) 

Comparing these results to stochastic simulations, see Fig.[5l 
we recover a good agreement. 

In the second case, the densities of both states are in 
pure LD or HD phases, i.e. boundary layers occur in both 
density profiles. As discussed at the end of Subsec. 14.21 
the latter implies aj, etj ^ 0, thus ff,£|. and £t>£| are 
found to be the localization lengths, describing how far the 
boundary layers extend into bulk. We are thus interested 
in their scaling behavior. 

Both d ] \d l are different from 0; approaching the de- 
localization transition line translates into d\ er — > 0. Our 
starting point is again the solution to the expanded Fokker- 
Planck equation (fTTjl given in App. El From Eq. flU]) we 
infer that £ ~ (det^4) _1 . Together with a, ~ d l and 
Eq. (fT4|). this implies for the localization lengths £j ~ 
(d 1 )^ 1 . For the scaling behavior in L, we observe aj ~ L^ 1 
and find ^ ~ L~ x . Together, the scaling behavior is given 

by 

As discussed in the beginning of this subsection, the 
same scaling behavior emerges in TASEP for the localiza- 
tion length £ of the boundary layer when approaching the 
coexistence line. 



5 Summary 

The appearance of a domain wall derealization is a par- 
ticular feature of the exclusion process with internal states 
introduced in Refs. 13,23]. Indeed, in the simplest driven 
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exclusion process, the TASEP [2|[T6] . such a delocaliza- 
tion of a localized domain wall does not emerge. There, 
boundary layers characterize the different phases. A delo- 
calized domain wall forms at the phase boundary between 
low- and high-density phases, while localized domain walls 
do not build up. Upon coupling to external reservoirs [25l 
I26j . localized domain walls may appear in bulk, separat- 
ing a low-density from a high-density region. Varying the 
system's parameters, their location changes continuously, 
and they may leave bulk through the left or the right 
boundary, causing continuous transitions to pure LD or 
HD phases. Domain wall delocalization, corresponding to 
a discontinuous transition, does not emerge. 

The exclusion processes with internal states exhibits, 
in addition to continuous phase transitions similar to the 
ones described above, discontinuous transitions. In par- 
ticular, in phase space, two regions with localized domain 
walls may be separated by such a discontinuous transi- 
tion; there, the domain wall position jumps from a certain 
location in bulk to another. This jump is connected to a 
delocalization of the domain wall: Approaching the discon- 
tinuous transition, the domain wall delocalizes and, upon 
crossing the transition, re-localizes at the other position. 

We have presented a detailed study of this delocaliza- 
tion phenomenon. First, we have investigated the system's 
behavior at the discontinuous transition. From mean-field 
considerations, delocalized domain walls have been iden- 
tified which perform coupled random walks, resulting in 
smoothened density profiles that lack sharp shocks. Sec- 
ond, we have investigated the delocalization of a localized 
domain wall upon approaching the discontinuous tran- 
sition. Using the domain wall picture, we have set up 
a quantitative description of the delocalization. Starting 
from coupled random walks of the domain walls, we have 
expanded the corresponding Master-equation in the Kra- 
mers-Moyal formalism, used the linear noise approxima- 
tion by van Kampen and obtained an analytically solvable 
Fokker-Planck equation. The latter has revealed power- 
law dependences of the domain wall width a on the dis- 
tance d to the delocalization transition and the system 
size L: a ~ \Ld\~ 1 ! 2 . These findings have been validated 
by stochastic simulations. 

As in similar driven diffusive systems where a local- 
ized domain wall appears [25,26 , the domain wall's width 
tends to zero with increasing system size, and is propor- 
tional to its inverse square root. However, in the present 
system, the domain wall can delocalize, namely upon ap- 
proaching the discontinuous transition. There, the width 
diverges, being proportional to the inverse square root of 
the distance to the transition. 

We believe that the above discussed domain wall de- 
localization represents a robust and generic phenomenon 
that may emerge in other driven diffusive system a well, 
such as weakly coupled antiparellel transport |35j or in 
models for intracellular transport on multiple lanes that 
take motors' internal states into account [361(37] . Its identi- 
fication in other non-equilibrium systems will shed further 
light on their universal phenomenology. Its observation in 
real systems, such as intracellular transport of kinesins 



moving on the parallel protofilaments of a microtubulus, 
constitutes a challenge for future research. 
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A The Fokker-Planck equation and its 
solution 

In this Appendix, we want to give more details concerning 
the technical parts of Sec. [4] The coefficients of the Fokker- 
Planck equation (|SJ) are derived, and the solution to the 
resulting matrix equation (fT3|) is given. 



A.l The coefficients of the Fokker-Planck equation 

The Fokker-Planck equation for the distribution of the 
domain wall positions reads 

d t P(y, t) = -di[ai^)P% t)] + ^BtjdidjPiy, t) . (18) 

The coefficients and Bij are given by Eqs. ((9]), we exem- 
plify their calculation for <2j. According to ([9]), the latter 
is connected to the changes , arising from the processes 
(i), (ii), (v) and (vi) (see Subsec. 14. 2p . The changes to- 
gether with the respective rates are read off from Tab. [1] 
We obtain 

a T (y) =L~ 1 Aj 1 [ - a c T ff (l - a\ s ) + /3 c T ff (l - /&) 
l 



n p T [i- P l ]dx~ n j P l [\-p ] ]dx 

Jo Jo 



(19) 



where we used the definitions (jT5|) for the distances 

to the delocalization transition. The other coefficients are 

obtained along the same lines; together, they read 

a T (y) =ITMfV, 

ai (y) =irM^V, 

B n (y) =L~ 2 A- 2 {aHl - cJ) +^(1 - + f2j tot 

B n (y) = -L- 2 A^Al 1 {t2 J t ° t 

+ n J [ p i(x)-p L (x)] 2 dx}, 

B u (y) =L- 2 Aj 2 {a l {l - a 1 ) + (3 1 {1 - P l ) + nj tot , 

[p^x) -p^(x)] 2 dx\ . (20) 
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For the expanded Fokker-Planck equation (fTTj) , the matrix 
elements A^ are given as the derivatives Aij = dja,i(y) | _ Q : 



n 



Tl 



IT 



L- 1 A^ 2 (d r A r )d^ 
L-'AfidrA^ 



II 



L^Af^A^ 



L- l A^ 
IT 1 AT 1 A 

L- 1 



T i 



(21) 



The key point of the analysis is the observation that 
in the determinant of A, 

detA = L- 2 {A^ 1 [z\^(a T Z\ T ) + Aj\diA r )]d 1 
+ Aj 1 [Aj 1 (d l A l ) + A- 1 (diA l )]dt 

'■[(d^A^AO - (d^id^dU^ , 

(22) 



the terms independent of d) , d^ have canceled. Thus, for 
d^,d^ — > 0, we encounter detA — > 0. We show in the 
following that this causes the domain walls to delocalize 
when d T , d^ — ► 0. 



A. 2 Solution of the Fokker-Planck equation 

For the solution of the expanded Fokker-Planck equation (fTTjl 
we already anticipated the form 

P(y, t) ~ exp [ - \{Z-%rfyi + f^V] ■ (23) 

To solve the resulting matrix equation 

SA T + AS + B = , (24) 

for E ', we write it in vector form: 



2A 



it 



2A 



% A n + A u A n 







2A 



,r 




+ \B U =0. (25) 



Note that B and S are symmetric. The above equation 
can be inverted to give 



S 



TT 



2(A n + An) detA' 



detA + A^ -2A n A u 



-Ail^n 
A 2 



2A n A u 

-2A TT A iT dct A + A 




-An A 



(26) 



For the scaling behavior, we observe S ~ det A. When 
the derealization transition is approached, as det A — > 
0, also £ — ► 0. Thus ZJ^ 1 , describing the widths of the 
domain walls, diverges; and the domain walls delocalize. 
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